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摘要 : 为 了 快速 获取 副 反 射 面 随 天 线 俯仰 角 变 化 的 位 置 ,并 通过 调整 机 构 对 副 面 位 置 进 
行 修正 来 提高 天 线 综 合 性 能 , 提出 一 种 虚拟 双 目 视觉 测量 方法 , 通过 高 速 相 机 与 平面 镜 组 成 
的 虚拟 双 目 视觉 系统 对 天 线 观 测 时 副 面 位 置 进 行 模拟 实验 , 获取 被 测 物 随 实验 平台 空间 移动 
产生 的 位 置 变 化 图 像 ， 并 对 图 像 进 行 预 处 理 、 去 品 、 边 缘 检 测 、 特 征 点 提取 等 操作 ,重建 了 
被 测 物 随 实验 平台 移动 而 在 三 维 方向 上 的 运动 情况 ,重建 结果 在 三 个 方向 上 误差 均 小 于 1%， 
表明 该 方法 具有 较 高 的 测量 精度 ， 能 够 较 好 的 对 副 面 位 置 情况 进行 重 构 。 
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1 引 


随 着 科学 技术 的 发 展 , 射电 望远镜 观测 技术 取得 了 长 足 的 进步 ， 对 天 文学 的 发 展 起 到 了 
极 大 的 推动 作用 。 与 此 同时 ， 大 型 双 反 射 面 天 线 的 口径 、 频 段 也 在 不 断 增 加 ， 对 天 线 的 指向 
精度 、 增 益 等 指标 要 求 越 来 越 高 趾 。 天 线 副 面 位 置 与 主 面 位 置 的 匹配 对 天 线 指向 精度 而 言 起 
着 至 关 重 要 的 作用 。 但 在 天 线 实际 运行 过 程 中 ， 其 副 反 射 面 的 位 置 会 受到 风 荷 、 重 力 、 温度 
等 因素 的 影响 而 发 生变 化 ， 导 致 主 副 面 位 置 匹配 失调 ， 使 天 线 的 综合 性 能 下 降 ， 对 天 线 的 效 
率 及 指向 精度 产生 影响 四。 因此 ， 如 何 快 速 、 实 时 、 准 确 地 重 构 天 线 副 面 位 置 ， 对 提高 天 线 
指向 精度 有 十 分 重要 的 意义 。 

为 获得 高 精度 天 线 副 面 位 置 ， 王 锦 清 等 人 推导 并 仿真 了 副 面 位 姿 和 全 息 口 径 面 相位 的 
解析 关系 ， 通 过 口径 面相 位 数据 ， 计 算 获 得 上 海天 马 65 米 望远镜 (以 下 简称 “TM65m”) 
副 面 位 姿 随 俯仰 角 改变 的 变化 情况 Bl。 许 谦 等 人 基于 应 变 传感器 和 模 态 受 加 原理 ， 分 别 通过 
采集 天 线 发 生变 形 后 副 反 射 面 撑 腿 、 副 反射 面 自 身 应 变 值 重 构 了 副 反射 面 的 位 姿 内 。Martin 
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Siiss 等 人 采用 摄影 测量 ， 对 天 线 副 面 随 俯 仰角 变化 而 产生 的 位 置 变化 进行 测量 ， 获 得 了 x、 
y、z 三 个 方向 上 的 位 置 变化 曲线 。 该 方法 获取 数据 量 较 大 ， 对 计算 体系 要 求 较 高 ， 不 能 够 
获取 天 线 副 面 的 实时 位 置 申 。Tonino Pisanu 等 人 和 江 永 琛 等 人 都 采用 PSD 法 测量 了 天 线 玉 
面 随 俯仰 角 变 化 而 在 x、y 两 个 方向 上 产生 的 位 置 变化 。 该 方法 可 以 获取 天 线 副 面 在 x、y 
两 个 方向 上 的 高 精度 位 置 ， 但 在 z 轴 方 向 上 的 位 置 需要 通过 其 他 方法 来 获得 [571。 侯 雷雨 等 
人 采用 激光 跟踪 仪 测量 了 副 面 调整 机 构 在 0" 和 45? 条 件 下 副 面 的 位 姿 情况 ,该 方法 具有 较 高 
的 精度 ， 在 对 天 线 副 面 调整 机 构 的 高 精度 标定 中 起 到 很 好 的 作用 ， 但 代价 太 大 加。 尹 航 等 人 
基于 惯性 测量 和 授 入 式 技术 , 解决 了 大 型 射电 望远镜 副 面 摆动 较 大 无 法 精确 定位 问题 。 但 该 
方法 使 用 惯性 测量 技术 ， 会 造成 误差 累积 ， 系 统 误差 越 来 越 大 ， 需 要 及 时 调整 四。 

本 文 将 采用 虚拟 双 目 测量 系统 , 对 天 线 副 面 在 空间 中 移动 过 程 中 的 三 维 位 置 进行 模拟 测 
量 , 得 出 目标 物体 在 空间 运动 的 三 维 运动 曲线 。 该 方法 只 需要 对 目标 物体 进行 一 次 拍摄 ， 就 
可 以 获得 副 面 在 该 位 置 的 三 维 信息 ,不 需要 很 复杂 的 程序 ， 同时 不 会 造成 误差 累积 。 运用 该 
方法 不 仅 能 节约 成 本 ， 且 能 够 解决 双 目 视觉 左右 相机 同步 难 的 问题 。 


Le, 


2 虚拟 双 目 视觉 原理 


交汇 式 双 目 视 觉 是 基于 同一 物体 在 两 图 像 中 的 视差 以 及 两 相机 之 间 的 平移 和 旋转 关系 
对 该 物体 空间 三 维 位 置 进行 求解 中。 其 原理 如 图 1 所 示 。 为 降低 成 本 ， 同 时 提高 拍摄 同步 
性 , 使 用 两 个 平面 镜 与 一 台 高 速 相 机 组 成 虚拟 双 目 视觉 系统 , 该 方法 与 交汇 式 双 目 视觉 原理 


基本 相同 ， 其 原理 如 图 2 所 示 。 
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滨 Camera 
1 交汇 式 双 目 视觉 原理 图 2 虚拟 双 目 视觉 原理 
Fig.1 The principle of convergent binocular vision Fig.2 The principle of virtual binocular vision 


如 图 2，P(x, y, z) 为 空间 中 任意 一 点 ， 会 分 别 在 左 镜面 和 右 镜面 上 成 像 ,高 速 相 机 对 左右 
镜面 上 成 的 像 进行 拍摄 ， 对 所 拍摄 图 像 进行 分 割 ， 会 得 到 两 幅 关 于 点 P(x, y, 2z) 的 像 ， 相 当 于 
有 两 台 虚 拟 相 机 对 点 P(x, y, z) 进 行 拍摄 。 其 分 析 方 法 与 真实 双 目 视觉 分 析 方 法 比较 起 来 基本 
相同 09。 


Ni 


Ow-XwYwZw 为 世界 坐标 系 ，O1-X1Y1Z1 为 左 相 机 摄像 机 坐标 系 ，O2-X2Y2Z2 为 右 相 机 摄 
像 机 坐标 系 。P; 点 为 空间 点 Pw 在 左 相 机 所 成 像 中 的 像 点 , Ps 点 为 空间 点 Pw 在 右 相 机 所 成 像 
中 的 像 点 。 
双 目 视觉 测量 中 ， 世 界 坐 标 系 与 摄像 机 坐标 系 的 关系 为 : 


Xw 
u 1 / dx 0 | 0 0 0 
R 人 yw 
Zc|v|= 0 l1/dy 7vo | | 0 f 0 0 
0 1 Zw 
1 0 0 1|10 0 1 0 | (1) 
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式 中 : (nv) 为 像 点 在 图 像 像素 坐标 系 中 的 坐标 ，〈Xw YwZw) 为 目标 点 在 世界 坐标 系 中 
的 坐标 。K 为 摄像 机 的 内 参 和 矩阵 ， 由 fdx、fdy、n、v 组 成 ， 由 摄像 机 自身 决定 。R 为 世界 

> 坐标 系 与 摄像 机 坐标 系 之 间 的 旋转 矩阵 ，T 为 世界 坐标 系 与 摄像 机 坐标 系 之 间 的 平移 矩阵 。 

十 R、T 称 为 摄像 机 的 外 参 和 矩阵 。 若 要 求 得 目标 点 在 三 维 空间 的 准确 位 置 ， 就 必须 求 得 摄像 机 

的 内 参 和 矩阵 和 摄像 机 的 外 参 和 矩阵 ， 这 个 过 程 称 为 摄像 机 的 标定 。 

若 摄 像 机 标定 已 完成 ， 且 像 点 P:、P2: 的 像素 坐标 均 已 获取 ， 那 么 由 式 〈1) 可 得 : 
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式 中 : (Cuull) 、 (hzv021) 为 像 点 Pl、P2: 的 齐 次 像素 坐标 ， (xw yw zw 1) 为 目标 点 
P 的 齐 次 世界 坐标 。 对 式 (2) 、(3) 进行 拆 分 并 消去 未 知 数 z1/、z 可 得 到 关于 xw、yw、zw 
的 四 个 线性 方程 组 成 的 方程 组 : 


(Um311— mM) Xxw + (W11321 — m121) pw + (U11M332 — 1131)Zw = 1M141 — U11M341 
(Vim311 — m211) Xxw t+ (Vi1m321 — m221) yw + (V11m322 一 1710231)Zw = 11241 — V11m241 (4) 
(W211312 一 11112)Xw + (W21322 一 11122) yw + (M2711332 一 111132) Zw = 1M142 一 U21M342 


(V21312 — 212 )xw t+ (V21322 一 11222) yw + (V211332 — 11232)Zw = 111242 — V21N 342 


对 式 〈4) 采用 最 小 二 乘法 求解 ， 即 可 求 得 目标 点 在 空间 内 的 准确 三 维 坐标 口 。 
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3 实验 及 图 像 处 理 


= 


测量 系统 如 图 3 所 示 ， 
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一 台 高 速 相机 ， 两 个 平面 镜 ， 一 个 棋盘 格 标定 板 ， 以 及 高 精度 


四 维 平台 和 一 台电 脑 组 成 。 通 过 电脑 控制 平台 在 X、Y、2Z 三 个 方向 移动 ， 高 速 相 机 对 被 测 


物 运动 图 像 进行 采集 并 通过 数据 线 传 输 到 电脑 ,并 由 电脑 对 图 像 进行 处 理 计算 。 实 验 获 取 图 


片 如 图 4 所 示 : 


图 3: 测量 系统 


Fig.3 The measuring System 


图 4: 实验 获取 


图 上 


Fig.4 The pictures obtained from the experiment 


在 获取 实验 图 像 后 ,我 们 对 图 像 进行 分 割 , 分 别 得 到 左右 虚拟 相机 的 图 像 ， 并 采用 张 正 


友 标 定 法 对 摄像 机 进行 标定 ， 标 定 结果 为 : 


8355.82314 0 0 
Kl = 0 8436.68818 0 
一 74.28131 373.97025 1 
8478.41734 0 0 
Kr = 0 8482.76571 o| 
479.50000 539.50000 1 
0.99001 一 0.03307 0.13705 
R= | 0.03479 0.99934 “一 0.01019 
一 0.13662 0.01486 0.99051 


T= |[-225.66790 -8.68607 12.21499| 


经 过 摄像 机 标定 得 到 摄像 机 内 外 参数 ， 就 可 以 对 目标 物体 进行 测试 。 在 拍摄 过 程 中 , 上 


于 设备 自身 影响 ， 以 及 拍摄 环境 中 噪声 、 光 照 等 因 
干扰 ， 使 图 像 质量 下 降 。 故 在 得 到 图 片 之 后 ， 需 要 
质量 ， 尽 可 能 的 接近 目标 物体 的 真实 情况 。 


素 的 影响 ,会 对 拍摄 图 片 的 成 像 质量 造成 
对 图 片 进 行 一 系列 的 处 理 ， 来 提高 图 像 的 


以 其 再 图 像 为 例 ， 如 图 
机 的 标定 参数 ， 对 图 
素 信息 ， 将 特征 点 以 外 的 像素 值 


三 


多 


5 中 a 


多] 


如 图 $ 中 b 图 像 所 示 。 
接 下 来 对 图 像 进行 滤波 。 
有 椒盐 噪声 与 高 斯 噪声 。 采 用 高 


的 图 像 如 图 5 中 ec 图 像 所 示 。 


对 于 滤波 处 理 过 后 的 图 片 进行 
边缘 检测 结果 


为 特征 点 的 边缘 。 


如 


(b) 


图 


(a) 初始 图 像 


像 处 理 


(a) The original image 


边缘 检测 ， 
图 5 中 dd 图像 所 示 。 


片 进行 畸变 修正 ， 并 且 对 修正 
设置 为 0。 这 


(b) Image preprocessing 


像 所 示 为 虚拟 相机 所 得 到 的 


过 滤波 来 去 除 由 于 


斯 滤波 能 够 在 很 大 程度 上 保留 图 像 的 边 


提取 像素 值 发 4 


个 过 程 称 之 为 图 片 的 预 处 理 ， 


(c) 图 像 滤 波 


(c) Image filtering 


图 5: 图 像 处 理 


Fig.S The image processing 


突变 的 位 置 ， 


图 片 。 首 先 ， 


通过 摄像 


过 的 图 像 进行 掩 膜 处 理 ,保留 特征 点 的 像 
处 理 过 


后 的 图 像 


所 构成 的 


日 摄 环 境 中 光线 对 图 片 造成 的 影响 , 主要 
力 缘 信息 0。 滤 波 过 后 


| 线 即 


(d) 边缘 检测 


(d) Edge detection 


并 以 左 虚拟 相机 


对 于 经 过 一 系列 图 像 处 理 过 后 的 图 像 ， 对 其 特征 点 中 心 坐 标 进行 提取 ， 
的 摄像 机 坐标 系 为 世界 坐标 ， 由 式 〈4) 可 得 特征 点 在 空间 的 三 维 坐标 。 

模拟 实验 分 为 两 组 进行 ， 第 一 组 模拟 实验 两 平面 镜 距 离 被 测 物 0.5m 左右 ， 
实验 两 平面 镜 距 离 被 测 物 2m 左右 。 


在 第 一 组 模拟 实验 过 程 


，X 方 向 以 Smm 为 一 


个 步 长 ，y 方 向 以 Imm 为 


方向 以 6mm 为 一 个 步 长 ， 特 征 点 在 空间 三 维 位 置 坐标 重 构 结 果 如 表 1 所 示 : 
表 1: 特征 点 三 维 坐标 重 构 结 果 
Tab.1 The result of 3D coordinate reconstruction of feature points 
1 2 3 4 5 6 是 
x 50.437 55.477 60.525 65.570 70.609 75.638 80.669 
y 23.282 24.285 25.289 26.291 27.298 28.303 29.309 
z 1076.258 1070.216 1064.247 1058.299 1052.358 1046.386 1040.359 
由 三 维 坐 标 绘制 x、y、z 方 向 运动 图 像 如 图 6 所 示 : 


第 二 组 模拟 


一 个 步 长 ， 过 


85.621 


30.314 


1034.395 


已 知 平台 在 x、y、z 方 向 上 
际 运动 距离 之 间 的 误差 ， 计 入 
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Tab.2 The error between the calculated result and the actual motion distance 


The distanee of Platform movement (mm) 


The distance of platform movement (mm) 


时 : 坟 台 


The distance of platform maovement (Mmm) 


图 6: 坐标 变化 图 


Fig.6 The coordinate change diagram 


的 移动 步 长 ， 


结 果 与 实际 运 


表 2: 计算 结果 与 实际 运动 距离 误差 


1 2 3 
0.8% 0.96% 0.9% 
0.3% 0.4% 0.2% 
0.7% 0.51% 0.86% 


S0 


OH 


通过 相 邻 两 点 坐标 之 差 可 求 出 计算 结果 与 实 
动 距离 之 间 的 误差 如 表 2 所 示 : 


4 5 6 7 
0.6% 0.58% 0.62% 0.96% 
0.7% 0.5% 0.6% 0.5% 
0.98% 0.46% 0.45% 0.6% 


， 在 第 二 组 实验 过 程 中 选择 不 同步 长 进行 测量 ，x 方向 以 3mm 


为 一 个 步 长 ，y 方向 以 1mm 为 一 个 步 长 ，z 方 向 以 7mm 为 一 个 步 长 ， 特 征 点 在 空 
置 坐标 重 构 结果 如 表 3 所 示 : 
表 3: 特征 点 三 维 坐标 重 构 结果 
Tab.4 The result of3D coordinate reconstruction of feature points 
1 记 3 4 5 6 7 
x 153.921 156.943 159.968 162.986 166.009 169.022 172.048 
y 47.927 46.920 45.915 44.911 43.902 42.899 41.893 
z 2689.469 2682.529 2675.566 2668.563 2661.520 2654.457 2647.436 
由 三 维 坐 标 绘制 x、y、z 方向 运动 图 像 如 图 7 所 示 : 
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Fig.7 The coordinate change diagram 
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40.883 
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本 组 模拟 实验 计算 结果 与 实际 运动 距离 误差 如 表 4 所 示 : 
表 4: 计算 结果 与 实际 运动 距离 误差 


Tab.$ The error between the calculated result and the actual motion distance 


1 2 3 4 3 6 7 
Xx 0.73% 0.83% 0.59% 1.0% 0.43% 0.84% 0.59% 
y 0.7% 0.49% 0.48% 0.86% 0.24% 0.62% 0.99% 
元 0.85% 0.53% 0.04% 0.61% 0.89% 0.29% 0.83% 


由 两 组 模拟 实验 结果 x、y、z 方向 坐标 变化 图 上 能 够 看 出 ， 计 算出 的 x、y、z 三 个 方向 
上 的 坐标 变化 与 实验 平台 在 三 维 方向 上 的 运动 有 较 好 的 线性 关系 。 通 过 计算 结果 与 实际 运动 
距离 之 间 的 误差 , 可 以 看 出 相 邻 两 点 之 间 的 距离 与 实际 平台 运动 距离 之 间 的 误差 均 小 于 1%， 
能 够 较 好 的 符合 真实 情况 。 且 当 两 平面 镜 与 被 测 物 距 离 为 0.5m 时 求 得 结果 误差 与 两 平面 镜 
同 被 测 物 距离 为 2m 时 求 得 误差 均 小 于 1%， 证 明 该 算法 有 较 好 的 稳定 性 。 


让 


4 虚拟 双 目 视觉 与 真实 双 目 视觉 的 比较 


为 评价 虚拟 双 目 视觉 精度 与 真实 双 目 视觉 精度 的 差距 , 评估 虚拟 双 目 视觉 的 可 行 性 , 采 
用 真实 双 目 视觉 对 被 测 物 三 维 位 置 进行 测量 。 两 相机 与 被 测 物 距离 为 2m 左右 ,实验 平台 在 
X、Y、Z 三 个 方向 分 别 以 5mm、lmm、6mm 为 步 长 进行 移动 ， 通 过 对 图 像 的 处 理 与 一 系列 
的 计算 ， 得 到 被 测 物 在 空间 中 的 三 维 位 置 如 表 5 所 示 : 
表 5: 真实 双 目 视觉 特征 点 三 维 坐标 重 构 结果 


Tab.$ The results of 3D coordinate reconstruction of feature points in real binocular vision 


-本 


1 2 3 4 5 6 全 8 
X 10.728 13:7353 20.777 25.800 30.820 35.841 40.863 45.882 
y 32.112 33.117 34.122 35.126 36.131 37.133 38.139 39.143 
Z 2044.055 2050.086 2056.115 2062.145 2068.172 2074.197 2080.221 2086.244 


计算 结果 与 实际 运动 距离 之 间 误 差 如 表 6 所 示 : 
表 6: 真实 双 目 视觉 计算 结果 与 实际 运动 距离 误差 


Tab.8 The error between the calculation result of real binocular vision and the actual motion 


distance 
1 2 3 4 5 6 7 
x 0.50% 0.48% 0.46% 0.40% 0.42% 0.44% 0.38% 
y 0.51% 0.49% 0.45% 0.47% 0.44% 0.43% 0.37% 


元 0.51% 0.48% 0.50% 0.45% 0.41% 0.40% 0.38% 


对 比 虚拟 双 目 视觉 计算 误差 与 真实 双 目 视 觉 计算 误差 ， 两 者 均 小 于 1%， 虚 拟 双 目 视 觉 
与 真实 双 目 视觉 之 间 的 误差 之 间 的 差距 仅 为 0.3% 左 右 ， 说 明 虚 拟 双 目 视觉 能 够 较 好 的 达到 
真实 双 目 视觉 的 精度 ， 同 时 能 够 说 明 该 方法 具有 较 好 的 可 行 性 。 


5 结论 


本 文通 过 对 虚拟 双 目 视觉 原理 及 构成 系统 的 研究 , 通过 高 速 相机 与 平面 镜 组 成 的 虚拟 双 
目 视觉 系统 对 天 线 副 面 随 俯仰 角 变 化 而 在 空间 中 产生 的 位 置 变化 进行 模拟 实验 , 重建 了 被 测 
物 随 实验 平台 移动 而 在 三 维 方 向 上 的 运动 情况 , 结果 能 够 较 好 的 符合 真实 情况 。 同时 对 虚拟 
双 目 视觉 与 真实 双 目 视觉 进行 对 比 , 表明 该 方法 具有 较 好 的 可 行 性 , 为 快速 重建 和 修正 天 线 
副 面 位 置 提供 了 一 个 有 效 的 方法 。 


6 应 用 设想 


本 研究 方法 将 针对 大 口径 望远镜 因 不 同 工 况 导致 副 面 位 置 的 变化 , 进而 造成 天 线 指向 精 
度 及 效率 下 降 的 问题 , 通过 对 副 面 位 姿 的 实时 测量 为 后 续 的 修正 提供 重要 支持 。 该 测量 系统 
将 靶 标 固定 于 天 线 副 面 上 ， 两 平面 镜 固定 于 天 线 主 面 ， 高 速 相 机 固定 于 天 线 副 面 撑 腿 位 置 。 
调试 两 平面 镜 之 间 夹 角 以 及 镜面 与 高 速 相 机 的 相对 位 置 ， 保 证 部 标 能 够 在 两 镜面 内 完全 呈 
现 , 同时 保证 高 速 相机 能 够 对 两 平面 镜 中 的 成 像 情况 进行 记录 。 通过 实验 验证 测量 精度 达到 
0.01mm 量 级 ， 由 于 系统 简易 、 数 据 处 理 方法 成 熟 使 得 系统 长 期 工作 可 靠 性 高 ， 能 够 较 好 的 
满足 天 线 整 体 技术 要 求 。 
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Abstract: In order to quickly obtain the position of the subreflector changing with the elevation 
angle of the antenna, and to improve the comprehensive performance of the antenna by adjusting 
the position of the subreflector, a virtual binocular vision measurement method is proposed. The 
position of the subreflector of the antenna ls simulated by the virtual binocular vision System 
composed of high-speed camera and plane mirror, and the position change image of the measured 
object moving with the space of the experimental platform is obtained. The image preprocessing, 
denoising, edge detection，feature point extraction and other operations are carried out to 
reconstruct the movement of the measured object in the three-dimensional direction with the 
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movement of the experimental platform, and the reconstruction results have errors of less than 1% 
in the three directions. The results show that this method has high measurement accuracy and can 
better reconstruct the position of the subreflector. 


Keywords: Antenna subreflector; Virtual binocular vision; Position; Reconstruction 
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